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ABSTRACT 

We report the results of a high spatial (parsec) resolution HCO + ( J = 1 — ¥ 0) 
and HCN (J = 1 — >■ 0) emission survey toward the giant molecular clouds of 
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the star formation regions N 105, N 113, N 159, and N44 in the Large Magellanic 
Cloud. The HCO + and HCN observations at 89.2 and 88.6 GHz, respectively, 
were conducted in the compact configuration of the Australia Telescope Compact 
Array. The emission is imaged into individual clumps with masses between 10 2 
and 10 4 M Q and radii of < 1 pc to ~ 2 pc. Many of the clumps are coincident with 
indicators of current massive star formation, indicating that many of the clumps 
are associated with deeply-embedded forming stars and star clusters. We find 
that massive YSO-bearing clumps tend to be larger (>1 pc), more massive (M 
> 10 3 M ), and have higher surface densities (~ 1 g cm -2 ), while clumps without 
signs of star formation are smaller (<1 pc), less massive (M < 10 3 M Q ), and have 
lower surface densities (~ 0.1 g cm -2 ). The dearth of massive (M > 10 3 M Q ) 
clumps not bearing massive YSOs suggests the onset of star formation occurs 
rapidly once the clump has attained physical properties favorable to massive 
star formation. Using a large sample of LMC massive YSO mid-IR spectra, 
we estimate that ~ 2/3 of the massive YSOs for which there are Spitzer mid- 
IR spectra are no longer located in molecular clumps; we estimate that these 
young stars/clusters have destroyed their natal clumps on a time scale of at least 
~ 3 x 10 5 yrs. 

Subject headings: galaxies: individual (LMC) -- infrared: stars — instrumenta- 
tion: spectrographs — Magellanic Clouds — stars: formation — stars: evolution 



1. Introduction 



Despite their dominant role in shaping galactic structures and stellar content, the current 
understanding of massive star formation (M > 8 M Q ) remains incomplete. There are a 
number of reasons for this, many related to observational difficulties. Massive stars form 
deeply embedded in the densest regions of molecular clouds with high extinctions, m aking; 



them difficult to observe during their early formation stages (e.g-.lLada &: Ladall2003l ). Due 



to the nature of the stellar initial mass function (IMF; ISarpeterlll955f ) , they are rare obj ects 



and coupled with their short formation times (~ 10 5 yrs; e.g., see IChurchwelll 120021 . and 



references therein), identifying observa tional targets is diffi cult. Furthermore, high-mass 



stars are seldom found in isolation (e.g.. |Preibisch et al.ll200fl ). meaning the interpretation of 



observations is complicated by the presence of other cluster members, the local environment, 
outflows, winds, and ionizing radiation. 

In recent year s a rough evolutionar y scenario for massive star formation has emerged (se e 
recent reviews by iBeuther et al.l 120071 : iMcKee fc Ostrikerl 120071 : IZinnecker fc Yorkei 120071 ) . 
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Stars are born within giant molecular clouds (GMCs), which have radii on the order of tens 
of parsecs and masses that can exceed 10 4 M . These GMCs are hierarchical in nature, 
with substructures that span a range of scales, and the literature contains a c onfused lexicon 



of term inology describing this substructure. Following the nomenclature of IWilliams et al. 



(120001 ) . the largest (approximately parsec-sized) over-dense regions with typical masses of 
10 2 — 10 3 M are defined as clumps. A core is a smaller, denser region of molecular gas 
that will form a single star or a small multiple system. The details of how massive star 
formation proceeds from these basic units within the GMC is unclear. A protostar forms 
from gravitational collapse wi thin the core, and it may grow in mass by accretion of the 
immediate core material (e.g.. iMcKee &: Tan! 120031 ) . Alternatively, the star may grow from 



the accretion of clump material initially unbound from the core (e.g., iBonnell et al.lll997l ). 
with the most massive members of a cluster forming near the center of the clump where the 
gravitational potential and resulting accretion rates are highest. The Kelvin- Helmholtz con- 
traction timescale of a high-mass star is smaller than the star's accretion timescale, meaning 
accretion can still be occurring after the star has begun hydrogen burning. The radiation of 
the star will ionize its surrounding to an increasingly larger scale until a combination of ra- 
diation pressure, stellar winds, and outflows have dissipated the dense surrounding material. 
In this paper, we observationally explore the properties of clumps as they evolve from dense 
star-forming structures to the less dense medium dispersed by the newly-formed stars. 

Observations of dust and molecules can trace the core and clump material being con- 
verted into stars and provide direct insight into the initial conditions of star formation. There 
have been numerous quantitative studies of the pr operties of these GMC substructures con- 
ducted in Milky Way star formation regions (e.g ., IStutzki et al.lll988t IWilliams et al.lll994j ; 
Elmegreen fc Falgarondl 19961 ; Kramer et al.lll998l ) . These gas and dust surveys have revealed 
the dense clumps to have densities of n Ha ~ 10 5 cm -3 , gas temperatures of 10 — 30 K, radii 
of 0.1-1 pc, and gas masses ranging from a few hundred to tens of thousdands of M that 
are found to be consistent with a power law dN/dMocM a , where a is typically found to be 
between —1 and —2. The mass spectrum extends to lower-mass overdensities essentially 
unbroken over many orders of magnitude. Kinematic information about the dense gas is ac- 
cessible through the rotational transitions of molecules with large dipole moments(e.g., CS, 
HCN, HCO + ), and the massive clumps are found to be gravitationally bound structures. 

Given the spatial resolution limitations, the properties of the dense gas in extragalac- 
tic star formation regions is less const rained than in the Gal axy. Using the Swedish-ESO 
15m Submillimeter Telescope (SEST), iJohansson et al.l (119941 ) report the detection of eight 
molecules ( 12 CO, CS, SO, CN, HCN, HNC, SCO+, and H 2 CO) towards the LMC's N159 



HII complex, and lChin et al.l (J19971 ). who extended the observations to the N 113, N44, and 
N214 complexes, detected 12 CO, 13 CO, CN, CS, HNC, HNC, HCO+, HC 2 , and C 3 H 2 . While 
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the gas was unresolved in SEST's large beam (~ 50"), the presence of dense gas is implied 
by the detection of the emission lines of molecules with high critical densities. With multi- 
ple transitions, it is pos sible to estimate the gas temperature and density, and from 12 CO 
and 13 CO observations, iMinamidani et al.l (J201ll ) finds a temperatur e of T = 15 — 200K 
and density of n Ha = 10 3 cm" 3 in 32 13 CO clumps in 6 LMC GMCs. JBolatto et all (boooh 
finds a dense gas emission component of T=15 K and nn 2 = 10 5 cm -3 from 12 CO and 13 CO 
observations of N159 by using a two-phase gas model that also included a T=100 K and 
n H2 = 10 2 cm~ 3 phase. 

Such low spatial resolution observations guided several targeted observations of single 
clumps in order to more robustly determine their properties. Determining the total mass 
contained within the cold, dense phase implied by the above molecular detections is only 
p ossible if the dens e gas is spatially resolved. Higher resolution interferrometric observations 
bv lOtt et al.l (120101 ) of one of the emission peaks in N 159 resulted in the first detection of 
NH 3 in the LM C; the two detected lines implied gas temperature of ~ 16 K and a total gas 
mass of 10 4 M Q . IWong et al.l (120061 ) resolved a clump in N 113 with targeted HCN and HCO + 
observations, determining a clump virial mass of ~ 10 4 . Here we report the results of a high 
resolution study of dense molecular clumps in several star formation regions in the Large 
Magellanic Cloud (LMC) using HCO+ ( J = 1 -»• 0) and HCN (J = 1 ->• 0) as high-density 
gas tracers. The observations, the first systematic interferometric mapping of dense gas in 
the LMC, were conducted with the Australia Telescope Compact Array (ATCA). 

Some of the problems of studying star formation in the Galaxy can be mitigated in 
the LMC. Galactic star formation is concentrated along the plane of the Galaxy, where 
dust obscuration is high; however, the LMC is located at high Galactic latitude with little 
intervening Galactic material. Milky Way studies are also plagued by distance uncertainties, 
while the LMC is located a t a known distance of 50 kpc (for a review of recent distance 
calculations, see lAlved 12004 ) . making luminosity determinations relatively robust. One of 
the primary observational advantages to studying the LMC is its proximity relative to other 
extragalactic systems, which allows single massive YSOs (or tightly bound systems) to be 
resolved even in the mid-IR. Taking advantage of the ability to observe the entire galaxy, 
a number of studies have sought to identify an d catalog all the massive and intermediate - 
mass young stellar objects (YSOs) in the LMC Jwhitnev et al.ll2008l : IcTruendl fc Chull2009h . 
Combined with other studies that identified sites of on-going star formation via their infrared 
emission, maser activity, and compact Ha emission an d centimeter emission, the literature 



now contains an extensive catal og of LMC YSOs (e.g., Ilndebetouw et al.ll2004t I Chen et al. 



20091 . bOlOuEllingsen et al.ll2010l ) to which we will positionally compare the dense gas revealed 
by our observations in order to explore the relationship between star formation and the dense 
gas. 
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We have conducted HCO + ( J = 1 — V 0) and HCN (J = 1 — >• 0) observations with the 
ATCA of four active star formation regions in the LMC - N 105, N 113, N 159, and N44. 



With critical densities of > 10 5 cm" 3 at T=20K flSchoier et all 12003 ). HCO+ and HCN were 
selected to probe the high volume densities of molecular clumps. We have analyzed the 
observations to determine the properties of the dense gas and its association with ongoing 
star formation in the LMC. Section 2 describes our observations and data reduction proce- 
dure. Section 3 details how discrete regions of molecular emission are identified, how their 
properties are determined, and the distributions of these clump properties. In Section 4 we 
summarize the star formation activity within the four complexes and explore its relationship 
with the imaged molecular material. Finally, a discussion of the feedback from star formation 
into the clumps is presented in Section 5; included is a scenario relating the properties of 
the clumps to the evolutionary state of the young stellar population, which we use to derive 
the timescales involved in the massive star formation process. 



2. Observations 

The regions chosen for observation in this study were selected for their known association 
with massive star formation. N 105, N 113, N 159, and N44 are well-studied H II complexes in 
the LMC, each containing one to several GMCs with signatures of both recent (Ha emission) 
and on-going (e.g., maser and bright IR emission) star formation. For each complex, we 
targeted the 12 CO emission peak of the region's GMC(s) with the aim of imaging the dense 
gas associated with current star formation in the cloud. Figures 1-5 sho w recent single-dish 

12 



CO observation s from the Magellanic Mopra Assessment (MAGMA; iHughes et al.ll20f0l ; 



Wong et al.ll201ll ). a survey of the molecular content of the Large and Small Magellanic 
Clouds using the Mopra Telescope along with the 50% sensitivity of our ATCA observations 
toward the CO emission peaks. Note that with a primary beam full width at half maximum 
(FWHM) of ~ 33", to cover large sections of the molecular clouds, mosaics were composed 
consisting of multiple ATCA pointings. 

The 3 mm HCO + and HCN data were collected during 12 separate sets of observations 
between Septermber 2006 and September 2009. At the time, the ATCA correlator could be 
configured to observe two frequencies simultaneously in dual polarization mode, allowing us 
to observe both HCO+ (J = 1 -> 0; 89.1885 GHz rest frequency) and HCN (J = 1 -> 0; 
88.6318 GHz rest frequency), redshifted to the appropriate observational frequency. All 
observations were conducted using the same array configuration (H75), which has a maximum 
antenna spacing of approximately 90 m, giving a synthesized beam FWHM of approximately 
6". Table 1 lists the sky position of each region, dates of the observations, number of pointings 



-6- 



in the mosaic, HCO + synthesized beam sizes, typical achieved sensitivities (similar for both 
frequencies), and frequency channel width of the observations. Note that we have separated 
N44 into two distinct, non-contiguous regions: N44's northern molecular cloud, which we 
refer to as N44 Region 1, and the southernmost molecular cloud, N44 Region 2. 

The data were calibrated and images create d using the Austr alia Telescope National 
Facility (ATNF) version of the MIRIAD package (ISault et al.lll995l ). Flux calibration was 
performed using Uranus, PKS 1921-293, and Mars for N105, N113, and N44, respectively; 
N159 was flux calibrated using PKS 0537-752, PKS 1253-055, and Mars over its 5-day 
observation. The overall flux calibration uncertainty is estimated at 20%. Quasar PKS 
0537-752 was used as the phase calibrator for all observations. PKS 1921-293 was used to 
calibrate the bandpass for N 105 and N113; PKS 1253-055 was used for N 159 and N44. 
After imaging the maps were CLEANed to a 2-cr level. These cleaned data cubes are then 
used in the clump- identification procedure described below. The contours and gray scale 
pixels in the bottom panels of Figures 1-5 are th moment (integrated intensity) maps of the 
HCO + and HCN data cubes, masked with a smooth version to reduce noise. 

To test the effect of the observations' uv coverage on the recovered source flux, we 
modeled an observation of a Gaussian source using the mosaic pattern and uv-coverage 
of the N105 mosaic. The model source was placed near the center of the mosaic, and fake 
visibilities were generated and imaged using the same procedure outlined above for the actual 
data. We varied the full-width-half-maximum (FWHM) of the model source to be between 
3" and 21"; > 80% of the flux is recovered for sources with FWHM< 9" and fluxes can be 
measured within a factor of 2 accuracy for intrinsic sizes of FWHM=14" or smaller. Larger 
sources tend to be resolved out. A poor signal to noise ratio reduces the recovered flux 
fraction roughly similarly for all size scales. By the definition of source radius defined in the 
following section, we are sensitive to structures of radius R< 2.6 pc at the distance of the 
LMC. 



3. Identification and Characterization of Dense Molecular Clumps 

The HCO + and HCN observations show that the molecular material within N 105, N 113, 
N 159, and N 44 is clumpy with emitting clouds of varied sizes, shapes, and velocity structures 
(Figures 1-5). HCO + and HCN had previously been det ected in N113's GMC with lower 



angu l ar resolution single d ish observations with th e SEST ([Johansson et al.lll994t IChin et al. 



19971 ; iHeikkila et al.lll999l ) and Mop ra telescopes ([Wong et al.ll2006l ). However, the implied 
density of the GMC is ~ 200 cm -3 ( IWong et al.ll2006l ). significantly lower than the critical 
densities of HCO+ (1 - 0) and HCN (1 - 0) (1.8 x 10 5 and 1.2 x 10 6 cm" 3 , respectively, at 
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20K; from the Leiden Atomic and Molecular Database. ISchoier et al.ll2005l ). The implication 
is that the molecular cloud is clumped into substructures, and our mosaic ATCA obser- 
vations directly reveal this dumpiness. We investigate the variations in clump properties 
quantitatively by determining the sizes, linewidths, and masses of the clumps. 



3.1. Identifying Clumps 

We use the HCO + ( J = 1 — > 0) data as the primary clump tracer. The HCN data has a 
lower overall signal-to-noise ratio (SNR) than that of HCO + because of the typically weaker 
HCN flux. To identify HCO + clumps in the ATCA data we use the algorithm described in 



Rosolowsky fc Leroyl (120061 ) by using the software version of the proc edure available as an 



I PL p ackage called CPROPS. The procedure is described in detail by iRosolowsky fc Leroy 



(120061 ) and in the CPROPS users guide, but a description of the basic steps of the algorithm 
follows. A mask is first generated from the data to isolate significant emission by estimating 
the noise variance (cr rms ) and identifying emission cores of four adjacent channels with flux 
density > 2.5<r rms . Each core is then expanded to include all connected emission above the 
2cr rms level. These masking levels (called the 'threshold' and 'edge' values in CPROPS) can 
be user-defined, and we explored a range of values. We find that threshold = 2.5a rms and 
edge = 2<7 rms were able to distinguish significant emission from noise while also identifying 
faint clumps. 

Once the images have been broken into these 'islands' of significant emission, CPROPS 
assigns the emission to individual clouds. Molecular clouds are known to be clumpy, and 
the islands identified here often contain substructure themselves. When a single surface 
brightness contour contains multiple emission peaks, CPROPS attempts to separate the 
island into multiple clumps. The peaks are separated into distinct clumps if each clump is 
larger than 1 resolution element (the synthesized beam), the contrast between each peak and 
edge of the clump exceeds a threshold, and the individual clumps have significantly different 
measured properties (e.g., mass, velocity width) compared to the combined island. When 
separating the islands into individual clumps with CPROPS, we chose to use the ECLUMP 
keyword that allows emission shared within a single brightness contour level by two separate 
clumps to be partitioned between the clumps. The ECLUMP keyword invokes a 'friends- 
of-friends' algorithm such that the pixels immediately surrounding a clump (friends) are 
incorporated into that particular clump, pixels immediately surrounding these friend pixels 
(friends of friends) are then incorporated, and so on, until all pixels are assigned to a clump. 
Pixels that are 'friends' with multiple clumps are assigned to the clump of nearest projected 
distance. CPROPS identified a total of 46 HCO + clumps in the five regions - five in N 105, 



six in N 113, twenty in N 159, nine in N44's Region 1, and six in N44's Region 2. 

Because of the signal-to-noise ratio (SNR) threshold invoked to isolate significant emis- 
sion, the identification of clumps is sensitive to the depth of the observations. Requiring a 
clump emission peak to have at least 4 consecutive channels above 2<r rms requires all clumps 
with broad velocity dispersions (er v S> channel width) to have a peak SNR of more than 
~ 3 within a single velocity channel. Clumps with narrow velocity widths (i.e., with ve- 
locity dispersions comparable to or smaller than the spectral resolution) require a higher 
peak SNR in order to have 4 consecutive channels above the 2<j rms threshold. For a clump 
with a narrow velocity FWHM of 1 km s _1 , observations with channel widths of 0.2 and 
0.4 km s _1 would require clump core detections of 3.2a rms and 3.9o" rms in a spectral channel, 
respectively. These detection thresholds amount to a slightly different minimum detectable 
flux limit for each region. The minimum detectable HCO + flux for an unresolved clump with 
a velocity FWHM of 1 km s' 1 is Fhco+ — 0-19 Jy km s _1 in N 105, whose observations have 
the lowest noise. On the other hand, N113 and N44 Region 2 have the highest minimum 
detectable flux of Fhco+ — 0.26 Jy km s _1 . These detection thresholds are very similar and 
lower than all save one clump's flux, and we therefore assert that despite slightly different 
noise levels and channel widths, the observations in each region are sensitive to the same 
clump population. 

While most of the HCO + emission peaks are strong detections with a high spectral 
SNR (up to SNR^ 10), those with peak SNR's of ~ 3 may be described as weak and/or 
possible non-detections. Although our threshold requirement of 4 consecutive channels above 
2cr rms generally does a good job of separating noise from true sources, we allow for the 
possibilty of false detections by labeling weakly-emitting peaks as 'candidate clumps;' their 
classification as real clumps will require deeper imaging. Candidate clumps, indicated by the 
"c" supercript in Table 2, are identified by narrow line widths (which could be mimicked by 
noise) and a low spectral SNR at clump center and across the entire clump (Figure 6). We 
assert all other detections are strong enough to indicate the emission is real and originates 
from true dense clumps. 



3.2. Derivation of Clump Physical Properties 

The radii of the clumps are determined by CPROPS, and in Table 2, we list the results. A 
clump's major and minor axes are first identified using principal component analysis (PCA; 
identifies the axis along which there is the most variation), and the second moments along 
each of these axes is then determined. The rms radius is the geometric mean of these two 
second moments. This measured size is necessarily an underestimate due to the sensitivity 
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limits of the observation; any flux belonging to the clump below the detection threshold will 
not be used in the above calculation. To account for this missing flux, CPROPS measures 
the size of the clump as a function of minimum surface brightness (essentially varying the 
observation's noise) and extrapolates the size of the clump to a value with zero noise. This 
method can also help mitigate the problem of varying sensitivity limits between different 
observations. The synthesized beam is deconvolved from the resulting rms size, which is 
then converted to clump radius by assuming a particular brightness profile and distance. If 
the extrapolated clump size from the cleaned image is smaller than the synthesized beam, 
it cannot be deconvolved. Furthermore, marginally-resolved clumps suffer from substantial 
radii uncertainties. As can be seen in Figures 1-5, many clumps are either unresolved or 
are comparable in size to the beam (~ 1 — 1.5 pc at the distance of the LMC). This is 
consistent with the sizes of dense molecular clumps in Galactic massive star forming regions 
with measured sizes of ~ 0.1 — 1 pc. Given the angular resolution limitations of extragalactic 
3 mm observations, at the distance of the LMC, only the largest (> 1 pc) clumps are resolved. 

A visual inspection of the ATCA data shows that the HCO + emission for each c lump 



emanates from a somewhat larger region than the HCN emission. IWong et al.l (120061 ) esti- 
mates the HCN-determined radius of N 113's clump 4 is 89% of that determined from HCO + 
and finds a decreasing HCO + /HCN ratio on longer interferometer baselines. As larger base- 
lines probe smaller scales, both findings suggests that the HCO + emission originates from a 
larger region than HCN. Such a radial dependence could be explained if HCN probes higher 
densities than HCO + . Indeed, HCN ( J = 1 — > 0) has a critic al density of w fri t = 1-2 x 10 6 at 



20 K, compared to n crit = 1.8 x 10 5 for HCO+ ( J = 1 -> 0) f lSchoier et al.ll2005f ). Moreover 



the HCO + abundance may be reduced in the high density centers of cl umps due to its par 



ticipation in proton transfer reactions with abundant neutral molecules ( jHeaton et al.lll993l ). 
and the abundance of HCO + may be enhanced in the outer, photon-dominated clump layers 
via chemical reac tions involving C + , which is expected to be abundant in regions exposed 



to UV radiation (jHeikkila et al.lll999l ). Each of the H II complexes studied here contain a 
number of UV-producing O and B-type stars, plausible sources of ionizing radiation that 
could enhance the HCO + abundance in the outer parts of the clump. In the latter case, the 
size determined from HCO+ in a high ionizing radiation field region is larger than would be 
found with less ionizing radiation. 

CPROPS calculates the velocity dispersion a v (the velocity size) of each clump from the 
second moment along the velocity axis of all clump pixels. Like the rms clump size, the 
velocity dispersion is extrapolated to what would be expected from a zero-noise observation 
and deconvolved with the channel width. For a Gaussian line profile, the full-width half- 
maximum (FWHM, Av) line width is then given by Av = a/8 ln(2) a v and is presented in 
Table 2. 
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Numerous previous studies of massive star forming clumps indicate that the masses of 



the clumps can be estimated via the virial theorem (e .g., lLarsonlll98ll ; ISolomon et al.lll987l ; 
Saito et al.l 12006c IWong et al.1 12006c iMuller et al.ll2010l ). Assuming the clumps are spherical, 
have radial density pr ofiles described by a power law, and are virialized, the virial mass is 
given by the formula (ISolomon et al.l 119871 ) 



M, 



125M, 



0' 



(5 - 2/3) 
(3-/3) 



Av 2 R, 



where (3 is the power law exponent of the density profile (density oc r -/3 ; r is the distance 
from clump center), R is the clump radius in pc, and Av is the FWHM of the velocity line 
width in km sec -1 . The virial masses pr ovided in Table 2 assum e the C PRO PS-determined 
radius and a density profile with (3 = 1 (Ivan der Tak et al.l 120001 ); the provided uncertainty 
assumes only observational effects (uncertainties in Av and R) and does not account for the 
uncertainty of /3. Clumps with extrapolated sizes smaller than the beam have no radius 
estimates, and therefore M vir cannot be determined. 

The masses of the clumps can also theoretically be determined directly from the HCO + 
line fluxes if the HCO + abundance and the optical depth of the line can be estimated, which 
requires the excitation temperature T ex ) to be known. Because we only have one transition 
of HCO + , it is not possible to determine T cx from our data, so in order to estimate the 
masses of the clumps from the HCO + flux, we assume a single valu e of T ex =30K for all 



clumps, typical of massive clumps observed in dust continuum (e.g., iFaundez et al.l 12004 ; 
Beltran et al.l 120061 ). We will investigate the effect of choosing a different value for T ex be- 
low. We estimate r, the clump's optical depth, from r = — ln[l - Tb/(T ex - Tb g )], where 
Tfe is the clump's t ypical brightness te mperature, taken to be 1/2 the clump's peak bright- 
ness temperature (IBarnes et al.1 120111 ); T ex is the excitation temperature; and T bg is the 
background temperature (2.72 K). We can then estimate the clumps' masses under LTE by 
calculating the col umn density and int egrating; over the surface area of the clump. Following 
the formulation of IBarnes et al.l ( 1201 ll ). 



M LTE = 4.0 x 10 3 M n t dV 



(2) 



where V, the line velocity, is in km s _1 and Q, a clump's projected surface area, is in pc 2 . 
Because we are assuming an excitation temperature and HCO + abundance (n(HCO + )/n(H2) 
= 10~ 9 ; the upper limit of recent core chemistry models and a typical value derived for 
massive clumps; see IBarnes et al.ll201ll . and references therein), Mlte is only an estimate. 
Resolved clumps have M L te uncertainties of typically ~ 40% derived from the uncertainty in 
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clump size and velocity dipersion. Because calculating Mlte requires T^ and Q, its estimation 
for unresolved clumps is more uncertain. For the unresolved clumps, Table 2 reports Mlte 
by assuming the clump is the size of the beam; the uncertainty of Mlte for unresolved clumps 
is ~ 200%, and because the clumps' brightness temperatures suffer from beam dilution, our 
estimate is likely an underestimate. 

The Mlte uncertainties in Table 2 take into account the uncertainty of Q, T&, and Av, 
but does not consider the uncertaity in r resulting from T ex . T ex is not measurable with our 
data and is potentially the largest source of uncertainty in our M L te estimate. Clumps that 
truly have a higher excitation (possibly such as those with embedded YSOs) have a lower 
Mlte, while in low excitation clumps (e.g., those without embedded YSOs), Mlte is higher. 
The optical depths (and therefore M L te) at the lower T cx of 10K are at least a factor of 3 
higher than at 30K (and higher for the brightest clumps) and a factor of 2 lower at the higher 
T ex of 50K. The analysis presented here is with a uniform T ex for uniformity and simplicity. 



3.3. Distribution of Properties 

Table 2 documents the derived properties of the clumps including position, Av, radius, 
M vir , Mlte, and total HCO + and HCN flux within the clump. Figure 6 shows the .Fhco+, 
-Fhcn, and M L te distributions for the entire clump sample. In each panel, the histogram 
displays a binned representation of the distribution whose number counts are indicated by the 
left axis, while the points are associated with the right axis and indicate each distribution's 
corresponding cumulative function, which is independent of bin choices. The vertical dotted 
line shows the location of the sensitivity limit in N 113 and N44 for a clump with size equal 
to that of the synthesized beam and Av = 1 km s _1 ; this corresponds to a clump with Mlte 
of ~ 25OM . 

The mass function of a population of clouds is typically described by a power law fit 
to the differential mass distribution, dN/dM oc M a , where a is the power law index of the 
distribution. N (M'>M), the integral of the mass function above mass M, can then be 
described by N(M' >M) oc M a+1 . The power law index of the mass function for Galactic 
cores (0.5 M < M < 20 M ) has been found to vary betwe en —1.6 to —2.8 w ith a typical 
va lue of ~ 2 in loc al star formatio n regions such as P erseus (jEnoch et al.l 120061 ) . Ophiuchus 



( lYoung et al.l 120061 ). and Serpens ( Enoch et al.l 120071 ). Recent early results from the Planck 
satellite found a Galactic core and clump (10 M^ < M < 10 3 M Q ) mass distribution consistent 
with a = —2 ( IPlanck Collaboration et al.ll201ll ). We find the clump mass distribution to be 
consistent with a = —2, but the uncertainty of Mlte prevents a more detailed fitting. 
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4. Star Formation Activity within the GMCs 
4.1. Identifying and Characterizing Star Formation Activity within the GMCs 



Recently, two separate, extensive surveys of LMC star formation were conducted us- 
ing the Spitzer Space Telescope with the aim of identifying the galaxy's YSO population. 
We use the term YSO to describe any young star (typical ages of several 10 5 — 10 6 years) 
that is still surrounded by enoug h circumstellar materi al to have a substa ntial infrared ex- 
cess. T ogether, these two surveys, IWhitney et al.l ( 120081 . hereafter W08) and lGruendl fc Chu 
( 120091 . hereafter GC09) identified 1000+ sources described as candidate YSOs (or tightly- 
bound, young multiple systems). We have used these two surveys as the primary catalogs 
from which we determine the locations of YSOs within the H II complexes studied, but we 
have also combed the literature for further YSO identifications. In particular, detailed IR 



studie s of N44 by IChen et al.l (120091 ) and of N 159 by Ijones et al.l (120051 ) and IChen et al. 



( I2010l ) have ident ified several YSOs not identified in GC09 or W08, and numerous compact 
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H II region (e.g.. Ilndebetouw et al.ll2004t iMartin-Hernandez et al.ll2005l) an d maser surveys 



Lazendic et al 



2002 



Ellingsen et al.ll2010l ) have been used 



(e.g., 

to strengthen YSO identifications. Figures 1-5 indicate the locations of sources identified as 
YSOs in the literature and Table 3 lists their positions and properties. The aforementioned 
means of detecting YSOs are either flux limited or restricted to phenomena only associ- 
ated with the most massive YSOs. Therefore only massive and intermediate-mass YSOs are 
marked in the figures. 



Based upon SED fits to radiation transfer models of YSO systems, GC09 and I Chen et al. 



( 120091 ) suggest that the 8.0 //m flux of a YSO may be a good proxy for the central source's 
total luminosity, and therefore, its mass. Table 3 lists each YSO's 8.0 /xm magnitude; we 
classify those with [8.0] < 8.0 as massive YSOs, while those with [8.0] > 8.0 are intermediate- 
mass. GC09 demonstrated that their candidate YSO catalog was nearly complete (> 99%) 
in the explored color-magnitude space for sources brighter than the 8.0 magnitude threshold; 
we assume all massive YSOs have been identified within the imaged regions. 

For the purposes of this paper, the evoluionary state of the complexes' YSO population 
is of particular interest. Man y of the YSOs, p articularly the most luminous, have been 
confirmed spectroscopically by ISeale et al.l ( 120091 . hereafter SL09) with mid-IR Spitzer IRS 
spectra. SL09 classified the spectra based upon the dominant spectral features and purports 
these spectral groups represent the evolutionary development of a YSO from one highly 
embedded in its natal molecular material to more advanced-stage YSOs with emission typical 
of photodissociation and compact H II regions. The strength of the silicate absopriton feature 
may be among the more evolutionarily diagnostic spectral features, and the feature's strength 
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along with each YSO's spectral group is reproduced in Table 2. Note that the silicate 
strengths were determined using principal component analysis (PCA) and should not be 
interpreted as the optical depth (nor is it necessarily proportional to the optical depth) of 
the feature as the spectra are first normalized to a common luminosity (see SL09). Instead, 
it should be interpreted as a dimensionless measure of the depth of the silicate absorption 
relative to the continuum that can range between -1 (little to no silicate absorption) to 1 
(strong silicate absorptio n). Moreove r , in the deta iled analyses of the star formation activity 



within N44 and N159, I Chen et al.l ( 120091 . 120101 ) supplemented Spitzer data with optical 



and near-IR images to construct SEDs from ~ 0.3 /zm to ~ 100 /im. They proposed a 
classification scheme for YSOs based upon their SED shapes. It is suggested that these 
Types are evolutionary (Types I, II, and III are arranged youngest to oldest) and reflect the 
amount of obscuring circumstellar mat erial around the forming star. For the YSOs in N44 
and N 159 that were assigned Types by lChen et al.l (J2009i . l2010l ). we have listed their Types 



in Table 3. 



4.2. Physical Association Between Clumps and Star Formation 

In this study we have identified a total of 36 dense molecular clumps (and 10 candidate 
clumps) from their HCO + emission. Even a by eye inspection of Figures 1-5 reveals a strong 
correlation between the location of the YSOs and the clumps' dense gas. The relationship 
between the YSOs and clumps can be quantified by comparing the proximity of the clumps to 
the YSOs and the clumps to a hypothetical field of YSOs randomly distributed throughout 
the region. An object randomly placed in region N 159 (the region we have most completely 
imaged) has only a 10% probability of being located within 12" of an emission peak of a 
clump (i.e., only 10% of N159's imaged surface area is within 12" of a clump's emission 
peak). However, a majority of the YSOs identified in N 159 are found to be < 10" from the 
nearest clump's center, indicating their positions are significantly different than a random 
distribution and are correlated with the positions of the clumps. 

As described above, the borders of the clumps are defined by a two-sigma detection 
threshold. Of the 36 clumps identified in our ATCA observations, 25 contain what have 
been identified as YSOs within their borders, and 5 of those have associated maser emission. 
Of the 25 clumps containing YSOs, 19 have at least 1 high-mass YSO ([8.0] < 8), while the 
6 others have only intermediate mass ([8.0] > 8.0) YSOs. Only 1 of the candidate clumps 
contains a YSO. 

In many cases, the YSO positions as determined from their IR emission is offset from 
the emission peak of the clump determined from the th moment (integrated intensity) map 
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of the HCO + emission but is still located within the contours of the clump. Recent modeling 
has shown t hat massive YSOs wil l form at the densest centers of their natal parsec-scale 



clouds (e.g., iKrumholz et al.ll2010l ). so the YSOs significantly offset from the peak of the 
clump may be more evolved than those more coincident with the core. The hypothesis that 
a YSO located near the center of a clump is evolutionary younger than those located near 
the clump edge is observationally supported by our finding that most YSOs (4/5) with maser 
emission - thought to be a signpost of early massive YSO evolution - are located at or near 
the emission peaks of the surrounding clump. 

If the offset between YSO position and clump center is evoluti onary, then it m ight be 



expected for clump properties to change as a function of this offset. Hi ll et al.l (120051 ) noted 
a number of Galactic massive star-bearing cores/clumps with star formation tracers that are 
offset from the peak of the cloud, but find no correlation between the clump's properties and 
the distance between the star formation tracer and the clump's emission peak. However, 
because the clumps they studied have sizes that range between ~ 0.1 and 2 pc, the distance 
between the clump center and the YSO may be influenced by the size of the clump as a 
whole. In other words, small YSO-bearing clumps will necessarily have YSOs close to their 
centers. To remove the dependence on clump size, we have determined a 'distance factor' 
that takes into account both the projected YSO-clump center distance and the size of the 
clump. The distance factor, S, is defined as 

s =cTe- (3) 

where C is the projected distance between the YSO and the nearest clump emission peak, 
and E is the projected distance between the YSO and the nearest clump edge. E is defined 
as positive for YSOs within the borders of a clump and negative for YSO outside a clump. 
Under this definition, YSOs at the emission peak of a clump have distance factors of 0, those 
located directly on the edge have distance factors of 1, and YSOs outside of clumps have 
distance factors of > 1. For clumps containing more than one emission peak, Table 3 reports 
the smallest value of S. 

The HCO + clumps detected in this survey comprise four distinct classes (N, IM, CM, 
and EM) and are defined as follows: 

1. N (No YSO) : clumps with no observed indications of recent intermediate or high- mass 
star formation. 

2. IM (Intermediate-mass YSO): clumps containing (5 < 1.0) at least one intermedi- 
ate mass YSO and no high-mass YSOs. 
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3. CM (Center massive YSO): clumps containing at least one massive YSO that is 
closer to the clump's center than its edge, i.e. 5 < 0.5. 

4. EM (Edge massive YSO): clumps with at least one massive YSO that is close to 
the clump's edge, i.e. 0.5 < 8 < 1.0 and no massive YSOs with 5 < 0.5. 

In the following section, we compare the properties of the clumps within each of these 
classes. 



4.3. Relationship Between Clump Properties and Star Formation Activity 

For the purposes of comparing the properties of the clumps in each group, we restrict 
ourselves to the most robustly-determined properties, and thus exclude candidate clumps 
from the following analysis unless otherwise noted. Histograms comparing the Mlte distri- 
butions of the different classes of clumps (Figure 8) show that the starless clumps are in 
general less massive than those associated with current star formation. The clumps without 
YSOs range in M L te between 2.6 x 10 2 and 2.3 x 10 3 M , with a mean mass of 1.5 x 10 3 M Q ; 
YSO-bearing clumps have a mass range of 8.0 x 10 2 — 3.3 x 10 4 M Q and an average mass of 
7.0 x 10 3 M . The difference in mass distributions is amplified when clumps of class IM are 
excluded - clumps with high-mass YSOs (CM and EM) have an average mass of 7.9 x 10 3 
M Q . We used the Student's t-test to statistically determine the significance of the differences 
in the mass distributions between the clump types; Table 4 reports the results. For each 
pair of clump types, we list the probability that the two populations do not have different 
means. Values below ~0.05 (5%) indicate significant differences in population distribution 
means. Note that the hypothesis that the clumps with and without YSOs come from the 
same population is rejected to a > 99% confidence level. 

The dearth of low-mass clumps (M < 2 x 10 3 M ) that are currently forming high- 
mass stars suggests that the mass of a clump is directly related to its massive star-forming 
capability. This is certainly true at the most basic level as a star could not form from clump 
material with a mass exceeding that of the natal clump. But we find that even the least 
massive clumps in this study have ample material to form a massive star, indicating there 
is an inherent property to the most massive clumps that makes them more favorable to 
forming massive stars. Recall that because we are only sensitive to the most massive YSO 
population members, the clumps we identify as being without YSO may not truly be free 
of star formation, and may be forming lower mass stars that our star formation tracers are 
insensitive to. 
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The finding of a lower mass for clumps without massive YSOs supports the findings of 



Hill et al.l ( 120051 ). who studied a large sample of Galactic, 0.1 — 2 pc-sized 'cores' (clumps 
under the terminology adopted here) in massive star formation regions. Cold, dense cores 
of dust were detected using 1.2-mm continuum observations, and a gas-to-dust ratio was 
assumed in order to estimate the cores' masses. They found that cores without maser or 
UCH II emission (their primary tracers of star formation) have a median mass of 2.8 x 10 2 
M Q , while cores with either masers or UCH II regions have a median mass of 1.0 x 10 3 . 
Moreover, within their sample, clumps without star formation are on average smaller in size 
than those with. Because many of our clumps are unresolved, the size distribution of the 
clumps in our sample is less secure. However, we also find that starless clumps tend to be 
smaller than those with YSOs, and that significantly more YSO-less clumps (~ 80%) cannot 
be deconvolved from the beam compared to YSO-bearing clumps (~ 20%). 

We can estimate clump surface densities of our clump sample using Mlte along with the 
clump surface area; YSO-less clumps have surface densities on the order of 0.1 g cm -2 , while 
the clumps with massive YSOs have an average surface density of ~ 1 g cm -2 . Differences 
in surface densities may reflect differences in volume densities - YSO-bearing clumps (class 
CM, EM, and I) are denser than those without YSOs (class N). The higher volume density 
of YSO-bearing clumps is supported by the relatively high HCN/HCO + flux ratio seen in 
clumps with YSOs. For the reasons discussed in Section 3.2, HCN likely traces higher 
densities than HCO + , and because clumps radially decrease in density, the two molecules 
likely trace different volumes of gas. Therefore a higher HCN/HCO + ratio is expected from 
structures containing high volume density gas. Clumps with the lowest densities may have 
little to no HCN emission at all. Indeed, while the HCN/HCO + ratio is low for group N 
clumps (0.19), the average HCN/HCO + flux ratio for all YSO-bearing clumps is 0.27, and 
is even higher for CM clumps (0.36). 

In a recent study of the role of radiat i on fee dback on the fragmentation and star for- 



mation within gas clouds, iKrumholz et al.l ( 120101 ) find that only clouds with a high surface 
density (E > 1 g cm -2 ) are able to form massive stars. In a radiation-hydrodynamic sim- 
ulation of a 100 M cloud with an initial surface density of E = 0.1 g cm -2 (typical of 
Galactic low-mass star forming regions), the cloud fragments strongly, forming a number of 
low-mass stars and no massive stars. Conversely, clouds with higher surface densities (E > 1 
g cm -2 ) fragment less, and put most their mass into massive stars, even with a total cloud 
mass of only 100 M&. The o bservations presented here are consistent with the conclusions 



of the IKrumholz et al.l (120101 ) simulation - the ability of a clump to form a massive star is 



determined by its surface density (or by extension volume density). 

In summary, we find that the clumps forming YSOs are physically larger, more massive, 
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and denser than clumps that are not. How exactly the starless clumps fit into the picture of 
massive star formation is unclear. One possibility is that they are the evolutionary precursors 
to the clumps actively forming stars. They may represent the earliest stages of massive star 
formation, before the formation of a YSO massive enough to have been detected in any LMC 
YSO surveys. In this case, their lower masses, sizes, and density could be attributable to 
their youth, and they may grow in size, mass, and density if they accumulate more material 
from the molecular cloud in which they reside. Alternatively, it may be that these clumps are 
never destined to form massive stars at all. Distinguishing the two scenarios is difficult, and 
it may be that both hypotheses are true. For example, YSO-less clump 6 of N 159 is more 
massive than several of the clumps currently harboring massive YSOs, suggesting it could 
eventually form a star. Conversely, given the dearth of low-mass clumps with YSOs, clumps 
with masses of only a few 10 2 M may never form massive stars at all and may only form 
low-mass stars. In either scenario, because clumps with masses > 2 x 10 3 M almost always 
contain massive YSOs, the onset of massive star formation within these massive clumps 
occurs quickly; massive clumps do not have a long YSO-less stage. 



5. Clump Dissipation and Destruction 

5.1. Evidence for Clump Dissipation 

Along with a correlation between a clump's properties and the presence or absence of 
a YSO, there appears to also be a correlation with the location of a YSO within the clump. 
Clumps with YSOs more centrally located generally have larger masses than clumps with 
YSOs closer to their edges (Figure 8). A similar correlation is also apparent in the HCN flux 
of the clumps (Figure 9); all CM class clumps have high HCN fluxes and higher HCN/HCO 4 " 
ratios. On the other hand, little to no HCN is detected towards clumps of class N, IM, or EM. 
Recall that we find that nearly all the maser-associated YSOs, possibly the youngest in our 
sample, are located very close to their clumps' peaks. This then implies that a massive YSO 
is first formed at an emission peak, where the clump densities were high and a core collapsed 
to form a star. YSOs located at the edges of a clump could represent more advanced stages 
of evolution where the YSOs have either moved away from the dense cores they formed in or 
have dissipated parts of their immediate surroundings, leaving only the parts of the clump 
far enough away or with high enough densities to have not been destroyed. 

Two alternative, non-evolutionary explanations are possible for the differences in mass 
distribution between clumps with YSOs at the centers and those with YSOs at the edges. 
First, higher mass clumps could be more centrally peaked while intermediate-mass clumps 
are more flocculent. YSOs would then form preferentially in the higher-mass clumps where 
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the bulk of the gas is located and the densities are highest. In less massive clumps, the mass 
distribution may be less peaked, and YSOs could form in a more distributed fashion. This 
scenario also accounts for the low HCN flux in intermediate-mass clumps - because it traces 
higher densities, HCN would only be detected towards the most centrally-peaked, densest 
clumps' centers. The second possibility is that because the most massive clumps also tend to 
be larger, a YSO has a greater probability of being located near the center of a more massive 
clump simply because the inner region of a massive clump has more surface area. This could 
result in more massive clumps being classified as CM than EM. While the largest clumps are 
of the CM class, the size distributions of CM and EM clumps do not differ significantly (see 
Figure 10). Our finding that masers are preferentially found at the centers of clumps may 
not necessarily argue against these scenarios. There may be a bias towards forming masers 
only in the most massive clumps or the densest regions of clump. Indeed, the existence of a 
maser requires high col umns that can shie ld the masing molecular material from the star's 



ionizing radiation (e.g., iGenzel et al.lll98ll ). 



There are 20 YSOs within our ATCA-imaged area (38%) that are not located within 
dense clumps. Or, said more precisely, are not currently located within dense clumps, as they 
likely earlier formed within one. YSOs not located within clumps therefore could represent 
the more advanced stages of YSO evolution after a YSO has dissipated its surroundings. 
The dearth of clumps with masses < 10 3 M harboring massive YSOs suggests that after 
the onset of clump disruption, destruction of the clump proceeds rapidly; a clump that is 
undergoing the process of dispersal spends very little time in a low-mass phase. The result is 
essentially a bimodal distribution of massive YSOs - those in massive (M > 10 3 M ) clumps 
and those not in clumps at all. 

We find that the YSOs located outside clumps tend to be visible at near-IR and even op- 
tical wavelengths, consistent with little circumstellar obscuration. Conversely, most sources 
in clumps are not visible at optical wavelengths, suggesting they are surrounded by a sig- 
nificant enough co lumn to extinct the short wavelength emission. Indeed, we find that the 



Chen et all (120091 ) Type classification scheme is very well correlated with the association of 
a YSO with dense molecular clumps, is a good indicator of a source's 'embeddedness,' and 
can be considered evolutionary in origin. Figure 11 shows the proport i on of YSOs within 



and outside of clumps that belong to each YSO Type from I Chen et al.1 (120091 . |2010j). Note 
that YSOs within clumps tend to be of Type I, I/II or II while YSOs outside of clumps are 
of Type II, II/III, or III. 

A YSO's propensity to be located within a clump is also well correlated with the presence 
of silicate absorption at 10 /zm in its infrared spectrum. While the silicate spectral feature 
could potentially originate from any silicate-bearing dust particles located between the YSO 
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and the observer, the bulk of this material is located in the YSO's immediate circumstellar 
environment. In time, the YSO dissipates its surroundings, removing circumstellar gas 
and dust, thus weakening the strength of the silicate absorption feat ure. There is little 



contri bution to the silicate feature's depth from foreground Galactic dust. ISchwering fc Israel 



( 119911 ) find that the extinction from Milky Way dust in the direction of the LMC is low; 
Ay = 0.2-0.5, corresponding to a 10 /im silicate feauture optical depth no of 0.01-0.03. The 
r 10 of the silicate absorption present in the YSOs' spectra measure between ~ 0.2 and 2.0, 
meaning Galactic dust accounts for at most 15% of the silicate feature's optical depth. 

Many of the YSOs in this study were observed in SL09, who measured the strength of the 
silicate feature relative to the continuum using PCA. Table 3 reports the PCA-determined 
silicate strength for the SL09 YSOs, and Figure 12 shows the close correspondence between 
silicate strength (circumstellar dust) and HCO + flux (circumstellar molecular gas). Figure 
13 shows the distribution of PCA silicate strengths for YSOs within (5 < 1.0) and outside 
(5 > 1) of a clump. The silicate feature is strongest for sources that are located centrally 
within a clump; all YSOs with 5 < 0.5 have PCA silicate strengths > 0. This is in contrast 
with the YSOs located outside of clumps: 4/9 of their members have silicate strengths < 0. 
Clearly the strength of the silicate feature is dependent on the presence of a circumstellar 
clump. 

The correlation is highly suggestive that the silicate absorption feature and HCO + 
emission originate from the same medium; the silicate absorption may be originating from 
dust on the same scale as the clumps. Unlike interstellar ices, silicates do not require cold, 
dense, radia tion shielded envir onments; they are detected even towards diffuse interstellar 



clouds (e.g., iBowey et al.lll998l ). The silicates contributing to the 10 /im feature may not be 
isolated to the YSO's immediate envelope or cold core, and may reside in the parsec-scale 
clump or even in the larger-scale, less dense molecular cloud. 



5.2. The Time Scale of Clump Disruption 

If the strength of the silicate feature is indeed an accurate indicator of the presence 
of a circumstellar clump around a YSO, given a complete sample of LMC YSO mid-IR 
spectra, we can estimate the number of YSOs located within clumps and therefore the 
relative time a YSO spends in clump-embedded and clump-free evolutionary stages. The 
SL09 IRS spectral catalog consists of nearly every YSO in the LMC that meets the criteria 
for being a candidate massive YSO in the GC09 catalog, and so we use it to statistically 
estimate the clump disruption time. In particular, because the average silicate strength of 
the entire SL09 massive YSO catalog is ~0, lower than the average value for YSOs located 
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within clumps (~ 0.2), it can be inferred that many YSOs from SL09 are not in clumps and 
that massive YSOs must spend a non-zero fraction of their lives not located within dense 
clumps. We estimate the number of YSOs in the SL09 sample located within clumps using 
the procedure outlined below. 

We fit the distribution of silicate strengths for YSOs observed in clumps as a Gaussian 
with an average value of 0.23 and a FWHM of 0.5 (bottom panel of Figure 13). YSOs not 
in clumps have silicate strengths less than ~ 0.3, and we therefore make the assumption 
that all YSOs in SL09 with silicate strengths greater than 0.3 are within clumps. Note that 
the SL09 silicate strength distribution (top panel of Figure 12) appears to have a larger 
right-hand wing than left-hand wing, suggesting that the entire distribution is composed 
of two separate distributions - one with a higher average silicate strength than the other 
(i.e., YSOs within and outside of clumps, respectively). To this end, we fit the entire SL09 
catalog of silicate strengths using a summation of two Gaussians: 1) the first (the right 
Gaussian curve in the top panel of Figure 13) is representative of the clump-embedded YSO 
population whose central position and width is defined by the aforementioned Gaussian fit 
to YSOs observed in clumps and whose height is determined by a fit to all SL09 YSOs with 
large (> 0.3) PCA silicate strengths, and 2) the second is representative of YSOs not in 
clumps, whose height, central wavelength, and width are allowed to vary as free parameters. 
The implied distribution of YSOs not in clumps is the left Gaussian curve shown in the top 
panel of Figure 13, and has a central value of —0.10 and a FWHM of 0.33. 

It is suggested in SL09 that sources with PCA silicate strengths of greater than 0.2 are 
embedded, while those with weaker PCA strengths are not. However, these results suggest 
that there are a number of YSOs embedded in clumps with silicate strengths less than 0.2. 
Therefore the total number of SL09 YSOs still embedded in their natal clumps is likely larger 
than is implied in SL09 (20%). The Gaussian fits to the entire SL09 catalog presented here 
suggests that 34% of all the massive YSOs in SL09 are within clumps; i.e., a massive YSO 
spends ~ | of its life embedded in a clump during the YSO evolutionary stages probed by 
SL09. 

GC09 (the catalog from which SL09 is based) identified YSOs via their Spitzer photom- 
etry, meaning SL09 does not have the spectra for YSOs outside of the GC09 photometric 
selection criteria. Specifically, GC09 misses the youngest, most embedded YSOs (too dim 
at Spitzer wavelengths) and the most evolved YSOs that have substantially dissipated their 
surroundings (too blue in color to meet GC09 color selection criteria). Figure 14, a graphical 
representation of the stages in massive star formation, shows the time period explored by 
SL09, tsL09- The highly embedded stage before Spitzer detection is identified as £a, while 
the more evolved stage immediately following its identification with Spitzer is t D . We divide 
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tsL09 into clump-embedded and clump-free stages, £b and tc, respectively; from the analysis 



we present above, we estimate that £b 



x t 



SL09 



and Sa- 



fe 



The total clump disruption time, the time between when a massive star forms in the 
clump and when the clump would no longer be detectable by our ATCA observations, is 
represented in the timeline as £sf = ^a + ^b, where £a is the time between the formation of 
a massive YSO and its detection by Spitzer. Our determination of £b sets a lower limit on 
tsF such that £sf > | x ^slo9 Y rs - During t SF , the clump is able to form cluster members, 
but after the dense clump is dispersed, massive star formation - possibly all star formation 
- within the clump will cease. 

Determining £sf is difficult as there have been to date no published surveys of the 
entire LMC targeting the most embedded YSOs that are not detectable at near- or mid-IR 
wavelengths, those in t^. Comparing this star formation stage to low-mass YSOs, these 
objects are comparable to the so-called 'Class 0' sour ces that are s t rong blackbody emitters 
in the far-IR. To estimate £a, we use the findings of ISewilo et al.l (120 lOl ). who conducted a 
search for Class YSOs in a 2° x 8° strip of the LMC ( approximately 25% of the LMC's 
total surface area) with the Herschel Space Observatory (IPilbratt et al.ll2010l ). The survey 
identified ~ 80 YSO candidates that were not identified in Spitzer studies in this subregion 
of the galaxy Extrapolating to include the entire LMC, we estimate there may be ~ 350 
YSOs that are too embedded to have bee n detected in the G C09 and W08 Spitzer surveys 
but can be identified with Herschel. The ISewilo et al.l (120101 ) survey is likely sensitive to a 



similar YSO stellar mass range as the GC09 and W08 Spitzer surveys (M star > 5 M ), so 
the two surveys are likely probing the same population of YSOs at two different stages of 
evolution. We therefore use the ratio of the number of YSOs detected by Spitzer (~ 1300) 
to that predicted to be detected by only Herschel (~ 350) to estimate £a ~ \ x ^slo9 an d 

^SF = ^A + ^B ~ Y2 X ^SL09- 

To estimate a numeric value for £slo9i we determine the total mass of the LMC's Spitzer- 
detected massive YSOs and assume a form ation rate for that population. A Salpeter IMF 
(dN/dM = M~ a ; a = 2.35; JSalpeterl Il955l ) predicts massive stars (M > 8 M ) account for 
~ 10% of the total mass of a stellar population. Therefore, the formation rate of massive 
stars, MSFR, is 10% the star formation rate of all stars, SFR, i.e. MSFR = 0.1 x SFR. 
The total mass of JV massive (8 M Q < M < 50 M ) stars is M(> 8 ) w 16iV M in a 
Salpeter IMF population, so if one assumes a constant star formation rate, it then follows 
that it takes time 
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to form N massive YSOs, where SFR is the average total galactic star formation rate in 
M Q yr -1 . Together, the GC09 and W08 studies have identified iV=300 distinct massive 
YSO candidates, where massive is defined as above, [8.0] < 8.0. We infer that the LMC's 
massive YSOs have a combined mass of 4.8 x 10 3 M Q . The LMC's global SFR has been 
estimated previously by using the in t egrate d H-a, UV, and far-IR fluxes as star formation 



tracers. For example, ICalzetti et all (120071 1 find SFRs of 0.05 M yr" 1 and 0.14 M vr -1 
from the galaxy's total far-IR and H-a flux, respectively, while W08 reports a SFRs of 
0.17 M yr -1 from the total UV flux. We adopt a typical literature SFR of 0.1 M yr -1 . 
Substituting these values of N and SFR into equation (4), we estimate £slo9 ~5x 10 5 yrs 
and tsF ~3x 10 5 yrs. Note that because many of the massive YSOs identified in GC09 
and W08 are likely multiple systems containing several massive stars, we may underestimate 
Mj^ >8 , and consequently £slo9 and £sf, by a factor of a few. 

Massive stars disrupt their surroundings with ioni zing and non-i o nizing photons, stel- 



lar and pre-stellar winds, and eventually supernovae. iTan fc McKed ( 1200 ll ) simulated the 
destruction from ionizing radiation, stellar winds, and radiation pressure of a 0.1 pc, 10 3 
M Q clump by a forming cluster. They find a destruction time of 0.3 — 2 Myrs, depending 
on the level of internal substructure of the clump, similar to the values we estimate here 
observationally. The shortest modeled destruction time, 3 x 10 5 yr was found for a uni- 
form, quiescent cloud containing 3 massive stars (M star > 8M ). Longer clump lifetimes are 
possible if the clump itself is highly clumpy; if most of the clump's mass is in dense cores 
with velocities set by virial equilibrium, the ionized gas that is destroying the clump can 
be contained through what they call "turbulent mass loading" for 1.5-2 Myrs. The clump 
destruction time we estimate from the LMC sample is comparable to or may exceed the 
destruction time of a uniform clump if we underestimate £slo9 for the reasons outlines above, 
suggesting that turbulent mass loading may indeed be an important process to confining the 
ionized gas and prolonging the life of the clump. 

Important to the theory of massive star formation is how this clump lifetime compares 
to the free-fall time of the clump, t ff = ^/3ir/(32Gp) = 1.38 x 10 6 (ntf/10 3 cm- 3 )- 1/2 yr. 
For Uh = 10 5 cm -3 , typical of Galactic star-forming clumps and approximately equal to the 
critical densities of the dense gas tracers HCO + and HCN, this yields tff = 1.4 x 10 5 yr. Our 
approximation of the clump destruction time presented above is admittedly crude (assumes 
an LMC SFR, requires an IMF extrapolation, etc.), however it gives a clump lifetime after 
the formation of a massive star of at least 2 times the free-fall time. Without a source of 
support against collapse, the clump should collapse in a free-fall time, implying a source of 
support against collapse such as turbulence or magnetic fields is required to allow the clump 
to last the several tff estimated here. 



-23- 



5.3. Clump Star Formation EfRciency 

The star formation efficiency (SFE) of the star-forming clumps, defined as the fraction 
of total mass converted into stars over the lifetime of a cloud, can be estimated by dividing 
the total mass of YSOs by the sum of the YSO and clump masses. The summed M L te of 
the massive YSO-bearing clumps is 1.5 x 10 5 M . As described in Section 6.2.2, by counting 
the number of massive YSOs and assuming a Salpeter IMF, we can estimate the total stellar 
mass within the YSO-bearing clumps, ~ 4.5 x 10 3 M©. Because most - possibly all - of 
these YSOs are unresolved clusters containing several massive stars, the total YSO mass 
is more realistically a factor of several higher. From these total clump and YSO masses 
we estimate the SFE of the clumps to be > 3%. The SFE of > 3% is on the high end of 



SFE's estimated for entire GMCs, typicall y ~ 1% (e.g . jDuerr et al.lll982l ) and higher than 



the global LMC SFE of 1% estimated by iFukui et al.l ( 119991 ). Our estimate is more crude 



than Galactic determinations as we are unable to individually count cluster members, but 
is markedly larger than that of LMC GMCs as a whole, an expected trend as SFE should 
increase as one probes higher densities. Our clump SFE estimates are more typical of those 
estimated for Gala ctic embedded clust ers such as Serpens, Rho Oph, NGC 1333, and NGC 



2071 (~ 10%) (see lLada fc Ladall2003l . and references therein). 



6. Summary and Conclusions 

We have conducted a survey of the dense molecular material in active star formation 
regions in the LMC using the ATCA. The observations, the first systematic interferometric 
mapping of dense gas in the LMC, were centered on the peaks of the regions' GMCs, where 
we identified a total of 46 dense molecular clumps with the CPROPS IDL package. Much 
of the molecular material is associated with signposts of on-going intermediate and massive 
star formation including maser emission, compact H II regions, and bright IR point sources. 
We have categorized the clumps by the extent and nature of their star formation content, 
and find correlations between the categories and clump properties such as mass and size. 
Our primary conclusions are summarized below: 

1. The gas within a GMC is clumpy with substructures that are directly revealed by the 
high volume density tracers HCO + and HCN. The ratio of the luminosity of the two 
transitions, i 7 HCN/-^HCo+ decreases with increasing clump mass, column density, and 
star formation activity, with HCN emission generally only being detected from the most 
massive and densest clumps. The HCN/HCO + ratio is ~ 1/3 in massive clumps (M 
> 10 3 M Q ), and on average a factor of two lower in their less massive counterparts. The 
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least massive clumps display no detectable HCN emission. Supported by the smaller 
physical scale from which HCN emission emanates, HCN likely is tracing higher volume 
densities than HCO + , and its observation is restricted to the higher density centers of 
only the most dense clumps, i.e., those most capable of forming stars. 

2. Determined from either the virial theorem or an integration over the surface area of the 
estimated column density, the clumps' masses span several orders of magnitude from 
a few 10 2 M to ~ 3 x 10 4 M Q . A fit to the integral of the mass function is consistent 
with a power law distribution with index a = —2. We determine the clumps to have 
radii between <1 pc and ~ 2 pc and surface densities of 0.1 - 1 g cm" 2 . 

3. The clumps are found to contain varying levels of current star formation, and we 
categorized them into four distinct classes: clumps with no on-going star formation (N), 
clumps containing only intermediate-mass YSOs (IM), clumps high high-mass YSOs 
located near their emission peaks (CM), and clumps with high- mass YSO located close 
to the clump edge (EM). 

4. Clumps with and without signs of recent or current star formation differ in their phys- 
ical properties. Massive YSO-bearing clumps tend to be larger (>1 pc), more massive 
(M > 10 3 M Q ), and have higher surface densities (~ 1 g cm -2 ), while clumps without 
signs of star formation are smaller (<1 pc), less massive M < M1O 3 M , and have lower 
surface densities (~ 0.1 g cm -2 ). The implication is that the ability of a clump to form 
massive stars is determined by its physical attributes. 

5. The dearth of massive (M > 10 3 M Q ) clumps without signs of massive star formation 
suggests the onset of star formation is rapid within a masive clump and occurs on a 
timescale significantly shorter than the clump's lifetime. 

6. Clumps with centrally- located massive YSOs (CM) generally have masses 2-3 times 
that of clumps with massive YSOs only on their edges (EM). We suggest that the 
difference may be evolutionary, with YSOs being born at the centers of clumps but 
becoming displaced from the clumps emission peak with time, while removing 1/2 to 
2/3 of the clump mass by the time it is located at the clump edge. The lack of clumps 
harboring massive YSOs with masses < 10 3 M suggests that after the onset of clump 
disruption, destruction of the clump proceeds rapidly; a clump that is undergoing the 
process of dispersal spends very little time in a low-mass phase. The result is essentially 
a bimodal distribution of massive YSOs - those in massive clumps and those not in 
clumps at all. 

7. There is a strong correlation between the intensity of HCO + and HCN emission and 
the strength of the 10 /zm silicate absorption feature seen in the mid-IR spectra of 
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YSOs. Objects located within dense clumps have the strongest silicate features while 
those outside of clumps show little to no silicate absorption. The close correspondence 
indicates the absorbing particles responsible for the silicate feature may reside in the 
parsec-scale structures defined by the clumps. 

Using the large sample of LMC massive YSO spectra presented in SL09, we estimate 
from the strength of the sample's silicate features that 34% of the YSOs in the catalog 
are embedded in clumps. By adopting a global SFR for the LMC of 0.1 M yr _1 , we 
estimate a massive YSO spends at least 5 x 10 5 yrs in the evolutionary stages present 
in the Spitzer YSO catalogs and derive a total clump destruction time after the onset 
of massive star formation of at least 3 x 10 5 yrs. 
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4.23 


± 0.36 


1.73 ± 


0.35 


EM 


6 


05:40:04.8 


-69:44:25 


2.5 ± 0.5 


0.8 ± 0.4 


9.3 ± 6.2 X10 2 


2.4 ± 2.4 


X 10 3 


1.63 


± 0.18 


0.60 ± 


0.17 


N 


7 


05:39:36.7 


-69:45:36 


7.4 ± 0.2 


2.5 ± 0.1 


2.6 ± 0.2 X10 4 


3.3 ± 0.3 


XlO 4 


16.54 ± 0.59 


5.80 ± 


0.57 


CM, mascr 


8 


05:40:02.3 


-69:45:00 


3.4 ± 0.6 




< 2.0 ± 0.4 x 10 3 


1.2 ± 2.4 


X 10 3 


1.02 


± 0.08 


0.34 ± 


0.08 


N 


9° 


05:39:41.5 


-69:46:56 


2.3 ± 1.1 




< 9.0 ± 2.0 X 10 2 


0.7 ± 1.4 


XlO 3 


0.61 


± 0.05 


0.15 ± 


0.05 


N 


10 


05:39:40.7 


-69:46:33 


3.2 ± 0.7 




< 1.7 ± 0.4 x 10 3 


0.8 ± 1.6 


X 10 3 


0.72 


± 0.05 


0.13 ± 


0.05 


I 


11 


05:39:58.6 


-69:45:03 


3.2 ± 1.2 




< 1.7 ± 0.4 x 10 3 


1.2 ± 2.4 


X 10 3 


0.65 ± 0.06 


0.07 ± 


0.06 


N 


12 


05:39:41.7 


-69:46:11 


3.9 ± 0.3 


0.5 ± 0.2 


1.4 ± 0.6 X10 3 


4.4 ± 3.6 


XlO 3 


3.65 ± 0.27 


1.84 ± 


0.26 


CM 


13° 


05:39:31.3 


-69:45:24 


1.7 ± 1.3 




< 4.9 ± 1.1 x 10 2 


4.2 ± 9.1 


XlO 2 


0.45 


± 0.04 


0.04 ± 


0.04 


N 


14'' 


05:40:05.2 


-69:44:48 


1.9 ± 1.7 




< 6.0 ± 1.3 X 10 2 


2.8 ± 6.1 


XlO 2 


0.31 


± 0.03 


0.08 ± 


0.03 


N 


15 c 


05:39:37.0 


-69:46:15 


2.7 ± 1.1 


0.4 ± 1.3 


5.5 ± 18.4 X10 2 


0.8 ± 5.1 


XlO 3 


0.58 


± 0.05 


0.17 ± 


0.05 


N 


16° 


05:39:34.4 


-69:45:10 


1.1 ± 1.4 




< 2.2 ± 0.5 X 10 2 


3.0 ± 7.1 


XlO 2 


0.42 


± 0.04 


0.09 ± 


0.04 


N 


17 


05:39:48.0 


-69:45:16 


1.6 ± 0.9 




< 4.1 ± 0.9 X 10 2 


0.7 ± 1.4 


X 10 3 


0.83 


± 0.06 


0.16 ± 


0.06 


N 


18 


05:39:40.3 


-69:45:40 


1.8 ± 1.5 




< 5.3 ± 1.2 x 10 2 


0.8 ± 1.7 


X 10 3 


0.86 


± 0.05 


0.20 ± 


0.05 


N 


19 


05:39:45.0 


-69:45:10 


2.1 ± 0.4 




< 7.6 ± 1.7 x 10 2 


1.9 ± 3.9 


X 10 3 


1.56 


± 0.11 


0.43 ± 


0.11 


N 


20 


05:39:42.8 


-69:45:06 


2.0 ± 0.8 


1.2 ± 0.4 


9.2 ± 7.5 X10 2 


1.2 ± 0.9 


X 10 3 


0.88 


± 0.08 


0.10 ± 


0.08 


N 


N 44 Region 1 


1 


05:22:03.4 


-67:57:44 


4.5 ± 0.3 


1.4 ± 0.2 


5.2 ± 1.1 xlO 3 


4.5 ± 1.3 


X 10 3 


3.10 


± 0.15 


0.69 ± 


0.18 


I 


2 


05:22:05.7 


-67:57:51 


4.9 ± 0.4 


1.5 ± 0.2 


6.7 ± 1.4 XlO 3 


3.7 ± 1.0 


X 10 3 


2.28 


± 0.12 


0.33 ± 


0.14 


I 


3 


05:22:02.1 


-67:57:50 


3.3 ± 0.3 


1.4 ± 0.2 


2.8 ± 0.6 XlO 3 


3.5 ± 1.1 


XlO 3 


2.98 


± 0.18 


0.49 ± 


0.21 


CM 


4 


05:22:12.8 


-67:58:32 


5.2 ± 0.4 


1.0 ± 0.2 


5.0 ± 1.3 XlO 3 


3.5 ± 1.4 


X 10 3 


2.46 


± 0.14 


0.85 ± 


0.16 


CM 


5 


05:22:08.5 


-67:58:06 


2.3 ± 0.6 




< 8.8 ± 2.0 X 10 2 


0.9 ± 1.8 


X 10 3 


0.68 


± 0.07 


0.01 ± 


0.08 


I 


(i 


05:22:12.3 


-67:58:14 


4.3 ± 0.4 


1.6 ± 0.2 


5.6 ± 1.2 XlO 3 


5.5 ± 1.5 


XlO 3 


3.03 


± 0.17 


0.41 ± 


0.20 


I 


7 


05:22:07.8 


-67:58:27 


4.2 ± 0.4 


1.4 ± 0.2 


4.7 ± 1.1 XlO 3 


4.7 ± 1.4 


XlO 3 


2.65 


± 0.14 


0.36 ± 


0.16 


EM 



CO 

to 



Table 2 — Continued 



Clump ID 
Number 


RA (.12000) 

[h:m:s] 


Dec (J2000) 
[d:/:«] 


Av 
[km sec - ] 


CPROPS Radius 
[PC] 


M vir 

[M ] 


M LTE a 

[M ] 


HCO+ flux 

[Jy km sec ] 


HCN flux 
[Jy km sec ] 


Class b 


8 C 
9° 


05:22:00.4 
05:22:04.3 


-67:57:38 
-67:58:18 


1.1 ± 1.2 
2.6 ± 0.9 






< 2.0 ±0.4 x 10 2 

< 1.1 ±0.2 x 10 3 


2.5 ± 5.8 X10 2 
0.8 ± 1.5 X10 3 


0.44 ± 0.04 
0.54 ± 0.05 


0.01 ± 0.05 
0.03 ± 0.06 


N 
N 












N44 


R,egion 2 











05:22:54.3 
05:22:51.8 
05:22:56.8 
05:22:55.2 
05:22:58.3 
05:22:53.9 



04:29 
04:16 
04:10 
04:07 
03:48 
04:32 



3.9 ± 1.2 
1.7 ± 1.0 
2.3 ± 0.6 
3.2 ± 0.4 
2.0 ± 1.2 
1.2 ± 1.5 



< 2.6 ±0.6 X 10 3 

< 5.0 ± 1.1 X 10 2 

< 8.8 ± 2.0 x 10 2 

< 1.7 ±0.4 x 10 3 

< 7.1 ± 1.6 X 10 2 

< 2.2 ±0.5 X 10 2 



1.2 ± 2.5 X10 J 
2.6 ± 5.5 X10 2 
0.9 ± 1.9 Xl0 3 

1.6 ± 3.3 X10 3 

4.7 ± 9.9 xlO 2 
2.2 ± 5.2 XlO 2 



0.76 ± 0.06 
0.45 ± 0.05 
0.85 ± 0.10 
1.39 ± 0.12 
0.55 ± 0.06 
0.24 ± 0.04 



0.04 ± 0.09 
0.06 ± 0.07 
0.07 ± 0.16 
0.04 ± 0.19 
0.13 ± 0.10 
0.01 ± 0.06 



N 

N 
EM 
EM 

I 

N 



a Note that the Mlxe uncertainties quoted arc statistical and do not include uncertainties in the absolute HCO flux due to flux calibration (20%) or HCO 
abundance. 

Clump classes: N — No YSO in clump; CM — Centralized massive YSO in clump; EM — Massive YSO only on the edge of the clump; IM — Intermediate mass YSO 
in clump, no massive YSO in clump. Clumps with previously-identified mascr activity arc also indicated. 

c Candidatc clumps (Sec Section 3.1). 



CO 
CO 



Table 3. YSO Properties 



YSO ID 


R.A (J2000) 


Doc (J2000) 


Clump YSO is 


[8.0] 


IRS 


PCA silicate 


Distance 


HCO+ Flux Density 


Notes 




[h:m:s] 


[d:/:«] 


located in 




Group a 


strength 


factor S 


[Jy beam km s )] 






N 105 


050949.80-685402.1 


05:09:49.80 


-68:54:02.1 


2 


7.80 ± 0.08 






0.15 


1.38 ± 0.15 




050950.11-685349.4 


05:09:50.11 


-68 


53:49.4 




9.35 ± 0.11 






2.87 


-0.03 ± 0.15 




050950.12-685426.9 


05:09:50.12 


-68 


54:26.9 




8.16 ± 0.09 






3.22 


0.06 ± 0.15 




050950.53-685305.5 


05:09:50.53 


-68 


53:05.5 


4 


5.24 ± 0.05 


P 


0.18 


0.21 


3.00 ± 0.15 




050952.26-685327.3 


05:09:52.26 


-68 


53:27.3 


1 


6.82 ± 0.07 


PE 


0.57 


0.31 


3.81 ± 0.15 


mascr 


050952.73-685300.7 


05:09:52.73 


-68 


53:00.7 


3 


5.86 ± 0.08 


PE 


0.06 


0.50 


1.40 ± 0.15 




050953.89-685336.7 


05:09:53.89 


-68 


53:36.7 


1 


7.85 ± 0.07 


P 


0.40 


0.68 


0.88 ± 0.15 




N 113 


051315.73-692135.9 


05:13:15.73 


-69:21:35.9 




6.97 ± 0.06 


PE 


0.06 


3.20 


0.58 ± 0.22 




051317.30-692236.7 


05 


13:17.30 


-69:22:36.7 


3 








0.49 


2.88 ± 0.22 




051317.54-692208.5 


05 


13:17.54 


-69:22:08.5 


fi 


9.45 ± 0.13 






0.71 


0.51 ± 0.22 




051317.69-692225.0 


05 


13:17.69 


-69:22:25.0 


4 


5.62 ± 0.06 


PE 


0.23 


0.40 


5.69 ± 0.22 


mascr 


051318.26-692135.5 


05 


13:18.26 


-69:21:35.5 




7.98 ± 0.07 


PE 


-0.25 


2.04 


0.60 ± 0.22 




051319.14-692151.0 


05 


13:19.14 


-69:21:51.0 


5 


6.78 ± 0.07 


PE 


-0.03 


0.61 


1.43 ± 0.22 




051320.75-692151.4 


05 


13:20.75 


-69:21:51.4 




8.75 ± 0.10 






1.83 


-0.02 ± 0.22 




051321.43-692241.5 


05 


13:21.43 


-69:22:41.5 


2 


5.67 ± 0.07 


PE 


0.15 


0.42 


2.79 ± 0.22 




051325.09-692245.1 


05 


13:25.09 


-69:22:45.1 


1 


5.53 ± 0.06 


P 


0.32 


0.14 


5.60 ± 0.22 


mascr 


N 159 


053929.21-694719.0 


05:39:29.21 


-69:47:19.0 


3 


7.20 ± 0.07 


P 


0.36 


0.68 


1.63 ± 0.25 


Type I; maser 


053935.99-694604.1 


05:39:35.99 


-69:46:04.1 




6.84 ± 0.06 


PE 


-0.08 


1.05 


0.08 ± 0.25 


Typo I 


053937.04-694536.7 


05:39:37.04 


-69:45:36.7 


7 


6.60 ± 0.07 


PE 


0.68 


0.37 


3.57 ± 0.25 


Typo I 


053937.53-694609.8 


05:39:37.53 


-69:46:09.8 


1 


5.82 ± 0.06 


PE 


0.04 


0.32 


1.95 ± 0.25 


Typo I/II 


053937.56-694525.3 


05:39:37.56 


-69:45:25.3 


7 


6.13 ± 0.08 


PE 


0.37 


0.36 


3.94 ± 0.25 


Typo I 


053940.78-694632.1 


05:39:40.78 


-69:46:32.1 


10 


8.02 ± 0.07 






0.64 


0.79 ± 0.25 


Typo II 


053941.89-694612.0 


05:39:41.89 


-69:46:12.0 


12 


5.93 ± 0.06 


PE 


0.38 


0.25 


2.98 ± 0.25 


Typo I/II 


053943.74-694540.3 


05:39:43.74 


-69:45:40.3 


18 


9.91 ± 0.05 






0.96 


0.07 ± 0.25 


Type III 


053945.18-694450.4 


05:39:45.18 


-69:44:50.4 




7.44 ± 0.06 


PE 


0.05 


3.44 


-0.31 ± 0.25 


Type III 


053947.68-694526.1 


05:39:47.68 


-69:45:26.1 




8.14 ± 0.03 






2.37 


0.07 ± 0.25 


Typo II 


053951.60-694510.5 


05:39:51.60 


-69:45:10.5 






PE 




4.92 


0.48 ± 0.25 




053952.60-694517.0 


05:39:52.60 


-69:45:17.0 




8.14 ± 0.03 






2.31 


-0.06 ± 0.25 


Typo III 


053959.34-694526.3 


05:39:59.34 


-69:45:26.3 




6.54 ± 0.06 


PE 


-0.37 


3.90 


-0.04 ± 0.25 


Type II 


054004.39-694437.6 


05:40:04.39 


-69:44:37.6 


5 


5.40 ± 0.06 


PE 


0.04 


0.63 


2.64 ± 0.25 


Type III 


054009.49-694453.5 


05:40:09.49 


-69:44:53.5 




7.31 ± 0.07 


PE 


0.29 


1.55 


0.38 ± 0.25 


Type I/II 


N 44 Region 1 


052156.97-675700.1 


05:21:56.97 


-67:57:00.1 




10.03 ± 0.10 


PE 


0.06 


8.13 


-0.40 ± 0.37 


Type III 


052159.6-675721.7 


05:21:59.6 


-67:57:21.7 




10.09 ± 0.24 






3.25 


-0.06 ± 0.37 


Type III 


052201.9-675732.5 


05::22:01.9 


-67:57:32.5 




9.07 ± 0.25 






2.28 


-0.02 ± 0.37 




052202.0-675758.2 


05:22:02.0 


-67:57:58.2 




7.77 ± 0.03 






1.29 


-0.13 ± 0.37 


Type II 


052202.11-675753.6 


05:22:02.11 


-67:57:53.6 


3 


7.68 ± 0.70 


PE 


0.29 


0.48 


0.43 ± 0.37 


Type II 


052203.30-675747.0 


05:22:03.30 


-67:57:47.0 


1 


8.36 ± 0.11 


PE 


0.15 


0.49 


1.26 ± 0.37 


Typo I/II 


052203.9-675743.7 


05 


:22:03.9 


-67:57:43.7 


1 


8.98 ± 0.10 






0.49 


1.24 ± 0.37 


Type II 



CO 



Table 3 — Continued 





YSO ID 


RA (J2000) 


Dec (J2000) 


CI 


ump YSO is 


[8.0] 




IRS 


PCA silicate 


Distance 


HCO+ Flux Dens 


ity 


Notes 




[h:m:s] 


[d:/://] 




located in 






Group a 


strength 


factor 5 


[Jy 


beam km s" 


- 1 )] 






052204.8-675744.6 


05:22:04.8 


-67:57:44.6 






8.76 ± 0.20 








1.15 




-0.13 ± 0.37 




Type 


II/III 


052205.2-675741.6 


05:22:05.2 


-67:57:41.6 






8.61 ± 0.23 








1.62 




-0.68 ± 0.37 








052205.3-675748.5 


05:22:05.3 


-67:57:48.5 




2 


8.13 ± 0.04 








0.58 




1.12 ± 0.37 




Type 


II 


052206.28-675659.1 


05:22:06.28 


-67:56:59.1 






8.40 ± 0.06 








6.66 




-0.19 ± 0.37 




Type 


II 


052207.27-675819.7 


05:22:07.27 


-67:58:19.7 






8.00 ± 0.06 




PE 


-0.13 


1.30 




0.18 ± 0.37 




Type 


III 


052207.32-675826.8 


05:22:07.32 


-67:58:26.8 




7 


8.30 ± 0.08 








0.45 




0.84 ± 0.37 




Type 


II 


052208.6-675805.5 


05:22:08.6 


-67:58:05.5 




5 


9.01 ± 0.09 








0.15 




0.80 ± 0.37 








052211.86-675818.1 


05:22:11.86 


-67:58:18.1 




6 


8.95 ± 0.08 








0.67 




0.04 ± 0.37 




Type 


I 


052212.24-675813.2 


05:22:12.24 


-67:58:13.2 




6 


8.97 ± 0.10 








0.12 




1.55 ± 0.37 




Type 


II/III 


052212.57-675832.3 


05:22:12.57 


-67:58:32.4 




4 


5.08 ± 0.05 




SE 


0.44 


0.16 




3.38 ± 0.37 




Type 


I 












N 44 Region 


2 
















052251.62-680436.6 


05:22:51.62 


-68:04:36.6 






8.60 ± 0.09 








2.10 




0.09 ± 0.22 








052255.12-680409.4 


05:22:55.12 


-68 


04:09.4 




4 


7.64 ± 0.07 




PE 


-0.03 


0.48 




0.54 ± 0.22 




Typo 


II 


052256.79-680406.8 


05:22:56.79 


-68 


04:06.8 




3 


7.77 ± 0.07 




PE 


-0.16 


0.68 




0.79 ± 0.22 




Type 


II 


052257.55-680414.1 


05:22:57.55 


-68 


04:14.1 




3 


8.60 ± 0.07 








0.97 




0.40 ± 0.22 




Type 


I/II 


052259.0-680346.3 


05:22:59.0 


-68 


03:46.3 




5 


10.46 ± 0.14 








0.95 




0.60 ± 0.22 




Type 


II 


a IRS Spectral Classes 


S — Silicate absorption only; SE — Silicate abs 
sive YSO Types for N 159 and N 44 taken from 


orption and fine 


■st 


ructurc 


lines; P - PAH < 


'mission: PE - 
respectively. 


- PAH emission and fine 


structu 


re lines 


When indicated, ma 


Chen et al.l 12009) 


andlCb 


en et al.l 120101) 





00 



Table 4. Student's t-test Analysis of Clump M L te Distributions 



CM+EM CM+EM+IM 



N 


0.01 


0.004 


0.01 


0.03 


0.03 


CM 




0.06 


0.08 


0.37 


0.17 


EM 






0.98 


0.18 


0.26 


IM 








0.21 


0.30 



Note. — The probability that the two clump populations 
(given by the row and column) do not have different means. 
Candidate clumps arc not included in this calcualtion. Clump 
classes arc defined in Section 5. 



CO 
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5:10:00.0 5G.0 52.0 4&.0 03:44.0 2:10:00.0 56.0 52.0 43.0 00:44.0 

Right ascension Right ascension 




56.0 52.0 48.0 09:44.0 5:10:00.0 

Right ascension 



5H.0 52.0 48.0 09:44.0 

Right ascension 



Fig. 1.— Spitzer 5.8 /im, MAGMA CO and ATCA HCO+ and HCN images of N 105. (Top 
left) Gray scale 5.8 /im image overlaid with MAGMA CO contours with the positions of 
YSOs marked with circles. Coutours are at 4xn (n=l,2,3,4). (Top right) Gray scale 5.8 
/im image with the 50% sensitivity and observational boundary of the ATCA HCO + mosaic 
in thick contours; clump boundaries are indicated with thin contours. (Bottom left) The 
Qth momen t (integrated intensity) of the ATCA HCO + masked version of the data cube (to 
reduce noise) in both gray scale and coutour. Contour levels are at 1o x 2 n (n = 0, 1, 2, 3, 
4...). Note that the lowest contour represents a 2-sigma detection and that each successive 
contour represents a factor of two increase in signal-to-noise. (Bottom right) Same as the 
bottom left panel, but for HCN. 



-38- 



30.O 



30.O 




30.0 26.0 5:13:20.0 15.0 10.0 S0.0 25.0 5:13:20.0 15.0 10.0 
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Fig. 2. — Same as Figure 1, but for N 113. 
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Fig. 3. — Same as Figure 1, but for N 159. 
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Fig. 4. — Same as Figure 1, but for N44 Region 1. 
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Fig. 5. — Same as Figure 1, but for N44 Region 2. 
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Fig. 6. — HCO + spectra of the 46 clumps identified by CPROPS. The dark solid line shows 
the spectrum at the central peak of the clump's emission, while the gray shaded region 
displays the average spectrum for all emission associated with the clump. Candidate clumps 
are marked with a superscript "c" (see text). 
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Fig. 7. — Distribution of F H co+ , -^hcn, and Mlte for the clump sample. The histogram 
displays the binned distribution, while the points show the integrated distribution. The 
left panel displays the distribution of -F H co+ (solid line histogram) and Fhcn (dashed line 
histogram) for both the clumps and candidate clumps. The vertical dashed line shows the 
location of the sensitivity limit for a Av = 1 km s _1 clump in N 113 and N44. In the center 
Mlte panel the solid line histogram is for the clump population and the dashed line histogram 
is for the candidate clumps. The thick solid line in the right panel for the integrated Mlte 
distribution is a power law of dN/dMoc M" 2 . 
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Fig. 8. — The distribution of Mlte f° r the clumps separated by clump groups. Top: dis- 
tribution for clumps (solid line with gray fill) and candidate clumps (dotten line with no 
fill) of group N ; Middle top: clumps with YSOs within their borders (groups IM, EM, and 
CM). Massive YSO clumps (EM and CM) are denoted by dark gray, intermediate-mass YSO 
clumps (IM) by a lighter gray. Middle bottom: clumps in group EM; Bottom: clumps in 
group CM. 
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Fig. 9. — The distribution of clump HCN fluxes separated by clump groups. Panels are as 
in Figure 7. 
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Fig. 10. — The distribution of clump surface areas separated by clump groups. Panels are 
as in Figure 7. 
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Fig. 11. — Proportion of YSOs from I Chen et al.l (120091 . l2010l ) of each SED type within (top 
panel) and outside of (bottom panel) clumps. YSO Types are indicated and differentiated 
by shade of gray, while the length of each Types's bar corresponds to the relative number of 
YSOs of that type. 
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Fig. 12. — PCA silicate strengths for the YSOs in our ATCA imaged area as a function of 
on-source HCO + flux density. Black filled circles are massive YSOs centrally located within 
a clump, gray filled circles are massive YSOs located on the edge of clumps, open circles 
are YSOs not located in clumps, and triangles are intermediate mass YSOs located within 
clumps. The line is a linear fit to the data to guide the eye. Spearman's rank correlation 
coefficient, R, is given in the upper left of the figure. 
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Fig. 13. — Distribution of the strength of the 10/xm silicate feature as determined by PCA 
for the YSOs contained within the observed areas of N105, N113, N 159, and N44 (middle 
and bottom panels) compared to the distribution for all massive YSOs for which SL09 has 
Spitzer IRS spectra (top panel). The middle panel shows the distribution for YSOs in the 
ATCA-images area not in clumps, while the bottom panel shows the same for YSOs located 
in a clump along with a Gaussian fit to the distribution. The thin curves in the top panel 
show the fit of this Gaussian to the SL09 catalog's high silicate strength wing (right Gaussian 
curve, representative of the distribution of clump-bound YSOs) and a Gaussian fit to the rest 
of the distribution (left Gaussian curve, representative of the distribution of YSOs outside 
of clumps). The dark curve is the summation of these two Gaussians. 
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Fig. 14. — A graphical representation of the stages of massive star forcmation within a 
clump. Time proceeds from left to right. 



